SSU-HEP-11/05 



Relativistic description of the double P-wave 
charmonium production in e+e^ annihilation 

A. P. Martynenko 

Samara State University, Pavlov Street 1, 44^011, Samara, Russia and 

Samara State Aerospace University named after S.P. Korolyov, 
^ ■ Moskovskoye Shosse 34, 44^086, Samara, Russia 



< 



X 



A.M. Trunin 

Samara State Aerospace University named after S.P. Korolyov, 
Moskovskoye Shosse 34, 44^086, Samara, Russia 



■ ^ ■• On the basis of perturbative QCD and the relativistic quark model we calculate 

Q^' relativistic and bound state corrections in the production processes of a pair of V- 

^. wave charmonium states. Relativistic factors in the production amplitude connected 

D i with the relative motion of heavy quarks and the transformation law of the bound 

' state wave function to the reference frame of the moving "P-wave mesons are taken 
into account. For the gluon and quark propagators entering the production vertex 

^ ' function we use a truncated expansion in the ratio of the relative quark momenta 

^ - to the center-of-mass energy ^/s up to the second order. Relativistic corrections to 

, the quark bound state wave functions in the rest frame are considered by means of 

CN I the Breit-like potential. It turns out that the examined effects change essentially the 

' nonrelativistic results of the cross section for the reaction e"*" + e~ — )• /ic + Xcj at the 

■ center-of-mass energy ^/s = 10.6 GeV. 
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I. INTRODUCTION 

The large value of the exclusive double charmonium production cross section measured 
at the Belle and BABAR experiments [l|, [2| reveals definite problems in the theoretical 
description of these processes [3|-[5| . Many theoretical efforts were made in order to improve 
the calculation of the production cross section e~^+e~ — )■ J j'^+rjc. They included the analysis 
of other production mechanisms for the state J/"^ + r]c jsl, 0| and the calculation of different 
corrections which could change essentially the initial nonrelativistic result [S-l^. Despite 



the evident successes achieved on the basis of nonrelativistic quantum chromodynamics 
(NRQCD), the light cone method, quark potential models for correcting the discrepancy 
between the theory and experiment, the double charmonium production in e~^e~ annihilation 
remains an interesting task. On the one hand, there are other production processes of 
the V- and P-wave charmonium states which can be investigated in the same way as the 
production of 5-wave states. Recently the Belle and BABAR collaborations discovered 
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new charmonium-like states in e~^e~ annihilation |19l . |20[. The nature of these numerous 
resonances remains unclear to the present. Some of them are considered as a V- and V- 
wave excitations in the system (cc). On the other hand, the variety of the used approaches 
and the model parameters in this problem raises the question about the comparison of 
the obtained results that will lead to a better understanding of the quark-gluon dynamics 
and different mechanisms of the charmonium production. Two sources of the changing of 
the nonrelativistic cross section for the double charmonium production are revealed to the 
present: the radiative corrections of order 0{as) and relative motion of c-quarks forming the 
bound states. An actual physical processes of the charmonium production require formation 
of hadronic particles in final states (bound states of a charm quark c and a charm anti- 
quark c), for which quantum chromodynamics can not provide high precision description. 
Further investigation of charmonia production can improve our understanding of heavy 
quark production and the formation of quark bound states. 

This work continues our study of the exclusive double charmonium production in e~^e~ 
annihilation in the case of a pure P-wave (cc) quarkonium on the basis of a relativistic quark 



model (RQM) [IJ, l2lH24l|. Note that the term RQM specifies the approach in which the 
systematic account of corrections connected with the relative motion of heavy quarks can 
be performed. The relativistic quark model provides the solution in many tasks of heavy 
quark physics. It uses a number of perturbative and nonperturbative parameters entering 
in the quark interaction operator. All observables can be expressed in terms of these param- 
eters. In this way we can check the predictions of any quark model and draw a conclusion 
about its successfulness. At the same time the existence of a large number of different quark 
models which are sometimes very complicated for the practical use put a question about the 
elaboration of the unified model containing generally accepted structural elements. Another 
approach to the heavy quark physics which does not contain the ambiguities of the quark 
models was formulated in [25]. As any other model of strong interactions of quarks and 
gluons the approach of NRQCD introduces in the theory a large number of matrix elements 
parameterizing nonperturbative dynamics of quarks. To a certain extent the microscopic 
picture of the quark-gluon interaction resident in quark models is changed by the global 
picture operating with the numerous nonperturbative matrix elements. The improved de- 
termination of color-singlet NRQCD matrix elements for S-wave charmonium is presented 
in j26| . Their study evidently shows that the account of relative order f ^ corrections signifi- 
cantly increases the values of the matrix elements of leading order in v. The correspondence 
between parameters of quark models and NRQCD which can be established, opens the way 
for better understanding of quark-gluon interactions at small distances. In this sense both 
approaches complement each other and could reveal new aspects of color dynamics of quarks 
and gluons. Thus, the aim of this study consists in the extension of relativistic approach to 
the quarkonium production from Refs.fl^, 21, 22[ on the processes e"*" + e~ — )■ /ic + Xcj and 
determination of the interrelationship with the predictions of NRQCD. 



II. GENERAL FORMALISM 

We investigate the quarkonium production in the lowest-order perturbative quantum 
chromodynamics. The usual color-singlet mechanism is considered as a basic one for the 
pair charmonium production. We analyze the reactions e"*" -|- e~ hc + Xcj, where the final 
state consists of a pair of P-wave {Xco, Xci, Xc2) and he charm mesons. The diagrams that 
give contributions to the amplitude of these processes in leading order of the QCD coupling 



FIG. 1: The production amplitude of a pair of "P- wave charmonium states in e~^e~ annihilation. Vh^ 
denotes the 7^- wave meson he and T^^cj denotes the V-wscve meson Xcj- The wavy line shows the 
virtual photon and the dashed line corresponds to the gluon. T is the production vertex function. 



constant are presented in Fig.l. Two other diagrams can be obtained by corresponding 
permutations. There are two stages of the production process. In the first stage, which is 
described by perturbative QCD, the virtual photon 7* produces four heavy c-quarks and 
c-antiquarks with the following four-momenta: 
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pi,2=2^±p, (P-P) = 0; gi,2 = -Q±g, (g-Q) = 0, 

Lp(0,p) 



(1) 



where P{Q) are the total four-momenta, p = Lp(0,p), q = Lp(0, q) are the relative four- 
momenta obtained from the rest frame four-momenta (0, p) and (0, q) by the Lorentz trans- 
formation to the system moving with the momenta P, Q. The momenta pi^2 of the heavy 
quark c and antiquark c are not on the mass shell: P12 = — p^ = — p^ 7^ m^. 

Relation (1) describes the symmetrical escape of the c-quark and c-antiquark from the mass 
shell. In the second nonperturbative stage, quark- antiquark pairs form the final mesons. 

Let consider the production amplitude of the "P-wave vector state and P-wave states 
Xcj {J = 0, 1, 2), which can be presented in the form [3, [H, 24|: 
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(iq 
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:Sp {^l^p, P)rf (p, g, P, Q)^lJq, Qhu + ^l,{q, Q)rf (p, q, P, QMSp, PH} 



where a superscript P indicates the P-wave meson, as(4m^) is the QCD coupling constant, 
a is the fine structure constant and Qc is the c-quark electric charge, ri^2 are the vertex 
functions defined below. The production processes + e~ he + Xcj contain the quark 
bound states. The transition of free quarks to the (cc) mesons is described by specific wave 
functions. The relativistic P-wave functions of the bound quarks \E''^ accounting for the 
transformation from the rest frame to the moving one with four momenta P, Q, are 
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X5;(Q,5,)(l + {)2) 
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2m{e{q) + m) 2m 



(4) 



where the hat is a notation for the contraction of the four vector with the Dirac matrices, 
vi = P/Mh^, V2 = Q/M^^j] e-p{Q, Sz) is the polarization vector of the spin-triplet state Xcj, 
e{p) = a/p^ + and m is the c-quark mass. The relativistic functions (3)-(4) and the vertex 
functions ri^2 do not contain the (5(p^ — M^/4-|-m^). More complicated factor including the 
bound state wave function in the rest frame presented in Eqs.(3) and (4) plays the role of 
the (5-function. This means that instead of the substitutions M^^ = 2e(p) and M^^j = 2e(q) 
in the production amplitude we carry out the integration over the quark relative momenta 
p and q. The amplitude (2) is projected onto a color singlet state by replacing Vi{0)uk{0) 
with a projection operator of the form Vi{0)uk{0) = 5ik/\/S- The relativistic wave functions 
in Eqs.(3), (4) are equal to the product of the wave functions in the rest frame and the 



spin projection operators that are accurate at all orders in |p|/m [IJ, |2J]. The expression 
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of the spin projector in a slightly different form has been derived primarily in [27| in the 
framework of NRQCD. Our derivation of relations (3), (4) accounts for the transformation 
law of the bound state wave functions from the rest frame to the moving one with four 
momenta P and Q. This transformation law was discussed in the Bethe-Salpeter approach 
in 28 1 and in the quasipotential method in [29^]. We use the last one and write the necessary 



transformation as follows: 

^1/2, pa/ TjW \ nl/2. 1^/3/ ovy 



^r(p) = Dr '""{Rtm"' "'^«p)^r(p), (5) 



where is the Wigner rotation, Lp is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D^^'^{R) is defined by 

\Dl/,\RY^) = S-\p^,2)S{P)S{p), (6) 

where the explicit form for the Lorentz transformation matrix of the four-spinor is 




We omit here the intermediate expressions giving rise to our final relations (2)- (4) jl4i |21 
The presence of the S{p- P) function allows to make the integration over relative energy p° if 
we write the initial production amplitude as a convolution of the truncated amplitude with 
two Bethe-Salpeter (BS) meson wave functions. In the rest frame of the bound state the 
condition p° = allows to eliminate the relative energy from the BS wave function. The 
BS wave function satisfies a two-body bound state equation which is very complicated and 
has no known solution. A way to deal with this problem is to find a soluble lowest-order 
equation containing the main physical properties of the exact equation and develop a per- 
turbation theory. For this purpose we continue to work in three-dimensional quasipotential 
approach. In this framework the double charmonium production amplitude (2) can be writ- 
ten initially as a product of the production vertex function Fi 2 projected onto the positive 
energy states by means of the Dirac bispinors (free quark wave functions) and a bound state 
quasipotential wave functions describing the P-wave mesons in the reference frames moving 
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with four momenta P, Q. Further transformations use the known transformation law for the 
bound state wave functions to the rest frame (5). The physical interpretation of the double 
charmonium production amplitude is the following: we have a complicated transition of two 
heavy c-quark and c-antiquark which were produced in e"^e~-annihilation outside the mass 
shell and their subsequent evolution firstly on the mass shell (free Dirac bispinors) and then 
to the quark bound states. In the spin projectors we have 7^ — just the same as 

in the vertex production functions Fi 2. We can not say exactly whether the charm quarks 
are on-shell or not in the spin projectors (3)- (4) because we should consider these structures 
as a transition form factors for the heavy quarks from the free states to the bound states. 
At leading order in ag the vertex functions Ff 2(p, P', Q, Q) can be written as 



rf(p,P;?,g) = 7M 



T^''{p,P;q,Q) = 7^ 
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j.D^'^ih), (8) 

jpD^^k,), (9) 



where the gluon momenta are ki = pi +qi, k^ = P2 + q2 and P = s = {P + Qy — (p_ 
P-, p+ are four momenta of the electron and positron. The dependence on the relative 
momenta of c-quarks is presented both in the gluon propagator Dn^{k) and quark propagator 
as well as in the relativistic wave functions (3), (4). Taking into account that the ratio of 
the relative quark momenta p and q to the energy ^/s is small, we expand the inverse 
denominators of quark and gluon propagators as follows: 
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(10) 

(11) 
(12) 



In the expansions (10)-(12) we keep terms of third order in relative momenta p and q. Sub- 
stituting (10)-(12), (3)-(4) in (2) we preserve relativistic factors entering the denominators 
of the relativistic wave functions (3)- (4), but in the numerator of the amplitude (2) we take 
into account corrections of third order in |p|/m and |q|/m. This provides the convergence 
of the resulting momentum integrals. Then the angular integrals are calculated using the 
following relations: 
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(14) 

where P^/? = {Qap ~ ^2a^^2^), Rv{q) is the radial momentum wave function of P-wave char- 
monium states, £ju((5, Lz) is the polarization vector in orbital space. The integrals in (13) 



dq, 
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and (14) are convergent due to the presence of relativistic factors. In this work we do not 
make expansions of all relativistic factors containing the relative momenta p, q as in our 
paper [ij]. This gives us an opportunity to calculate the corrections of the second order in 
|p|/m, |q|/m working with the convergent integrals. Both approaches of the expansion can 
be used. In the first one we expand all factors, immediately obtain the divergent integrals 
and should find additional arguments to fix its numerical value jTsj. On the second way, 
which we take in this work, we have no the divergent integrals in the corrections of second 
order. Undecomposed relativistic factors in (13)-(14) can be considered as a natural cutoff 
of the momentum integrals. For a specific P-wave state, summing over Sz and in the 
amplitude (2) can be further simplified as js^] 



J2 (1, L,; 1, 5,1 J, J.)4JQ, L,)4^(Q, Sz) 



'al3 
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(15) 



where {1, Lz]l, Sz\J, Jz) are the Clebsch-Gordon coefficients. Calculating the trace in the 
amplitude (2) by means of expressions (3)-(4), (8)-(9) and the system FORM [Sli, we find 
that the tensor parts of four amplitudes describing the production of P-wave charmonium 
states have the following structure: 



SoAhc + Xco) = A^ti^a^v^v^el"^, (16) 
SiAhc + Xci) = Bivip{vi ■ e*^J{v2 ■ 4 J + B2V2pivi ■ e*^J{v2 ■ 4 J + B^e*^^^p{v2 ■ 4 J + (17) 

Xci ' ^hc) + BQV2i3{£hc ' ^Xcl)' 



S2,fi{hc + Xc2) 
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Cieap(3^v"v^v^{v2 ■ elj + C2eapp^v'^V2el^^ + C3£<,p/37^/i!^i ^2 + (18) 



where the coefficients Aj, Ci can be presented as sums of terms containing the factors 
u = M^^JiMh^ + M^J, K = m/{MH^ + M^J and C,, = c\p)cJ{q) = [{m - e{p))/{m + 
e{p))Y[{m — e{q))/{m + e(g))]-', preserving terms with i + j < 2, and = {Mh^ + M^^jY /s. 
Exact analytical expressions for these coefficients are sufficiently lengthy (compare with the 
results written in Appendix A of our previous paper j2l|), so we present them in Appendix 
A of this work only in approximate numerical form using the observed meson masses and 
the c-quark mass m = 1.55 GeV. 

Introducing the scattering angle 6 between the electron momentum pe and the momentum 
P of the he meson, we can calculate the differential cross section da/d cos 6 and then the total 
cross section a as a function of r^. We find it useful to present the charmonium production 
cross sections in the following form {k = 0, 1, 2 corresponds to Xco, Xa and Xc2)- 
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where the functions F, 
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(1) 



(k = 



0,1,2) are written exphcitly in Appendix B. The factors 
(^(2) = 1/16 are introduced for the convenience. 



(20) 



The parameters Ui can be expressed in terms of momentum integrals J„ for the states he 
and Xcj as follows: 
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'4, LOT = OOiLOi. 



On the one side, in the potential quark model the relativistic corrections, connected 
with the relative motion of heavy c-quarks, enter the production amplitude (2) and the 
cross section (19) through the different relativistic factors. They are determined in the final 
expression (19) by the specific parameters Ui. The momentum integrals which determine the 
parameters Ui are convergent and we calculate them numerically, using the wave functions 
obtained by the numerical solution of the Schrodinger equation. The exact form of the wave 
functions \E'q''(p) and \E'o''"^(q) is important for improving the accuracy of the calculation of 
the relativistic effects. It is sufficient to note that the double charmonium production cross 
section a{s) in the nonrelativistic approximation contains the factor \R'f^^{0)\'^\R'^^^{0)\'^ . 
Small changes of the numerical values of the bound state wave functions at the origin lead to 
substantial changes of the final results. In the framework of NRQCD this problem is clo sely 



related to the determination of the color-singlet matrix elements for the charmonium [25 



Thus, on the other side, there are relativistic corrections to the bound state wave functions 
\E'o''(p); ^o''''(q)- III order to take them into account, we suppose that the dynamics of a 
cc-pair is determined by the QCD generalization of the standard Breit Hamiltonian in the 
center-of-mass reference frame 132-34 : 



H = Ho + AUi + AU2, Ho = 2 Jp2 + m2 - 2m 



Ar + B, 



AUUr) 
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(23) 

(24) 
(25) 



where L = [r x p], S = Si + S2, rif is the number of flavors, Ca = 3 and Cp = 4/3 are 
the color factors of the SU(3) color group, 7^; ~ 0.577216 is the Euler constant. To describe 
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TABLE I: Numerical values of the relativistic parameters (20), (21), (22) in the double charmonium 
production cross section (19). 



Meson (cc) 






j^exp^ GeV 


i?;,(0), GeV5/2 


ioi or 


UJ2 or cjs 


XcO 




0++ 


3.415 


0.33 


-0.28 


0.13 


Xcl 


l3Pi 


1++ 


3.511 


0.20 


-0.18 


0.07 


Xc2 




2++ 


3.556 


0.13 


-0.08 


0.01 


he 




1+- 


3.525 


0.17 


-0.14 


0.04 



the hyperfine splittings in P-wave charmonium we add to the standard Breit potential the 
scalar-exchange and vector-exchange confining potentials obtained in 35N37l|: 



AKtft(r) = fv 
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- (26) 



where we take the parameter fv = 0.7 for optimal agreement with the experiment. For the 
dependence of the QCD coupling constant as(/i^) on the renormalization point /i^ in the 



pure Coulomb term in (23) we use the three-loop result [38 



An IQnbilnL 



+ 



QAtt 



6?(ln'L-lnL- 1) + 62 , L = ln(/xVA2), (27) 



whereas in other terms of the Hamiltonians (24) and (25) we take the leading order approxi- 
mation. The typical momentum transfer scale in a quarkonium is of order of the quark mass, 
so we set the renormalization scale fi = m and A = 0.168 GeV, which gives = 0.314 for 
the charmonium states. The coefficients bi are written explicitly in 38|. The parameters 
of the linear potential A = 0.18 GeV^ and B = —0.16 GeV have the usual values of quark 
models. Starting with the Hamiltonian (23) we construct the effective potential model based 
on the ^^hrodinger equation and find its numerical solutions in the case of P-wave charmo- 
The details of the used model are presented in Appendix C. Then we calculate 



mum 



39 



the matrix elements entering in the expressions for the parameters Ui (22) and obtain the 
value of the production cross sections at ^/s = 10.6 GeV. Basic parameters which determine 
our numerical results are collected in Table I. The comparison of the obtained results with 
the previous calculations 0,0, 40, 41 1 and experimental data [l|, [2| is presented in Table II. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this paper we have investigated the role of relativistic effects in the production pro- 
cesses of P-wave mesons (cc) in the quark model. At the calculation of the production 
amplitude (2) we keep relativistic corrections of two types. The first type is determined by 
several functions depending on the relative quark momenta p and q arising from the gluon 
propagator, the quark propagator and the relativistic meson wave functions. The second 
type of corrections originates from the perturbative and nonperturbative treatment of the 
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quark- ant iquark interaction operator which leads to the different wave functions ^l/o''(p) 
'^^''■'{q) for the P-wave charmonium states. In addition, we systematically accounted for 
the bound state corrections working with the observed masses of "P-wave mesons {Xcj, he). 
The calculated masses of P-wave charmonium states agree well with experimental values 
(42I I (see Table III). Note that the basic parameters of the model are kept fixed from the 
previous calculations of the meson mass spectra and decay widths 0, 0, 43, 44 1- The 
strong coupling constant entering the production amplitude (2) is taken to be = 0.24 in 
accordance with the leading order QCD relation at = 2m. 

Numerical results and their comparison with the previous calculation in NRQCD are 
presented in Table II. We have included in it also new numerical results (several numerical 
mistakes contained in 22|] were corrected) obtained on the basis of quark model (23)- (26) for 
the production cross sections of a pair of S- and P-wave charmonium states. The exclusive 
double charmonium production cross section presented in the form (19) is convenient for a 
comparison with the results of NRQCD. Indeed, in the nonrelativistic limit, when u = 1/2, 
K = 1/4, Ui = {i > 1), = 16m^/s, the cross section (19) coincides with the calculation in 
[3|. In this limit the functions Fg'^''(r^) transform into corresponding functions from j^. 
When we take into account bound state corrections working with observed meson masses, 
we get u = M^^J{Mh^ + M^^J ^ 1/2, k = m/{Mh, + M^^J ^ 1/4. This leads to the mod- 
ification of the general factor in (19) and the form of the functions Fq^"* in comparison with 
the nonrelativistic theory (see [3|). It follows from the numerical values of the parameters 
Wj, presented in Table I, that the relativistic corrections could amount to 10 -r- 30% in the 
production amplitude. In fact their influence on the value of the production cross sections 
become considerably larger in the case of reactions e"*" + e~ — )■ /ic + Xci,c2- Only due to rela- 
tivistic contributions to the production amplitude the cross section cr(e+ + e~ — )■ /ic + Xci) 
increases in two times and a{e^ + e~ he + Xc2) in four times in comparison with the 
nonrelativistic calculation. Opposite influence on the value of the cross sections is deter- 
mined by relativistic corrections to the bound state wave functions in the rest frame. Indeed, 
relativistic effects change considerably the values of the nonrelativistic parameters B^{0), 
which transform into -Rp(O) (20). Different values of the mass of c-quark and nonpertur- 
bative parameters -Rp(O) make difficult the direct comparison of our numerical results with 
predictions of NRQCD. Note that nonrelativistic results obtained in our quark model are 
the following: cr(xco + ^c) = 0.101 fb, (j{xci + hc) = 0.417 fb, cr(xc2 + ^c) = 0.026 fb (compare 
with predictions of NRQCD in fourth column of Table II). Nevertheless, we can state that 
in all considered reactions e"*" + e~ — )■ /ic + Xcj the account of all relativistic effects leads 
to the decrease of the nonrelativistic cross section obtained in our model. It is necessary to 
point out once again that the essential effect on the value of the production cross sections 
he + Xcj belongs to the parameters -Rp(O) (20), a^, m. Small changes in their values can 
lead to significant changes in the production cross sections. In our model the nonrelativistic 
value i?p(0) = 0.24 GeV^/^. Accounting for the potentials (23)-(26) which give the good 
mass splitting for P-wave charmonium states, we observe simultaneously the decreasing and 
splitting in the parameter B!p{0) (see Table I). As a result the nonrelativistic cross sec- 
tions cr(e+ + — 7- /ic + Xci) and cr{e~^ + e~ he + Xci) decrease in three and six times 
correspondingly and a{e^ + e~ ^ he + Xco) reduces approximately on 25%. 

We presented a systematic treatment of relativistic effects in the P-wave double char- 
monium production in e~^e~ annihilation. We separated two different types of relativistic 
contributions to the production amplitudes. The first type includes the relativistic v/c 
corrections to the wave functions and their relativistic transformations. The second type 
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includes the relativistic p / ^/s corrections appearing from the expansion of the quark and 
gluon propagators. The latter corrections were taken into account up to the second order. 
It is important to note that the expansion parameter p/ \fs is very small. In our analysis of 
the production amplitudes we correctly take into account relativistic contributions of order 
0{v'^/c^) for the P-wave mesons. Therefore the first basic theoretical uncertainty of our 
calculation is connected with the omitted terms of order 0(p^/m^). Since the calculation 
of the masses of "P-wave charmonium states is sufficiently accurate in our model (the error 
is less then 1 %), we suppose that the uncertainty in the cross section calculation due to 
the omitted relativistic corrections of order 0(p^/m^) in the quark interaction operator (the 
Breit Hamiltonian) is also very small. Taking into account that the average value of the 
heavy quark velocity squared in the charmonium is (f^) = 0.3, we expect that relativistic 
corrections of order 0(p^/m^) should not exceed 30% of the obtained relativistic contribu- 
tion. Strictly speaking in the quasipotential approach we can not find precisely the bound 
state wave functions in the region of the relativistic momenta p > m which gives near 30% of 
the total value a (19). Using indirect arguments related with the mass spectrum calculation 
we estimate in 10% the uncertainty in the wave function determination. Larger value of the 
error will lead to the essential discrepancy between the experiment and theory in the calcu- 
lation of the charmonium mass spectrum. Then the corresponding error in the cross section 
(19) is not exceeding 15%. The significant improvement in the calculation of the relativistic 



corrections to the double charmonium cross section a{e^e — )■ J/'^ + 'r]c) was obtained in [45 



in the nonrelativistic QCD factorization formalism. The essential refinement was connected 
with many factors including the resummation of a class of relativistic corrections and the 
contribution that arises from the interference between the relativistic corrections and the 
corrections of the next to leading order in ag- In our work the appearance of divergent inte- 
grals over p and q for the corrections of order 0(p^/m^) and 0(q^/m^) is the consequence 
of expansions (10)-(12) used by us in order to perform analytically the angular integration 
in (2). The omitted corrections can be included and the obtained results can be improved if 
we calculate all integrals over the relative momenta p and q in (2) without any expansions. 
Another important part of the total theoretical error is related with radiative corrections 
of order which were omitted in our analysis. Our approach to the calculation of the 
amplitude of the double charmonium production can be extended beyond the leading order 
in the strong coupling constant. Then the vertex functions in (2) will have more complicate 
structure including the integration over the loop momenta. Our calculation of the cross 
sections accounts for effectively only some part of one loop corrections by means of the Breit 
Hamiltonian. So, we assume that the radiative corrections of order 0{as) can cause 20% 
modification of the production cross sections. We have neglected the terms in the cross sec- 
tion (19) containing the product of J„ with summary index > 2 because their contribution 
has been found negligibly small. There are no another comparable uncertainties related to 
the other parameters of the model, since their values were fixed from our previous consid- 



eration of meson and baryon properties j2J, |43|. Our total theoretical errors are written 



explicitly in Table II. To obtain this estimate we add the above mentioned uncertainties in 
quadrature. 
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TABLE II: Comparison of the obtained results with previous theoretical predictions and experi- 
mental data. 



State 

ZJ ZJ 


O'BABAR'X 
^rH2^charged>2 

(fb) [2] 


orH2^charged>2 

(fb) 


(^NRQCD 

ffU\ [Ql 

(tbj \6\ 


a (fb) 
[/ii 


a (fb) 
[40| 


a (fb) 

I A T 1 
[41] 


Our result 
(fb) 


J/"^ + XcO 


10.3 ±2.5^1:^ 


6.4± 1.7± 1.0 


2.40 ± 1.02 


6.7 


14.4 


17.9(6.35) 


14.47 ±5.64 


J/^ + Xcl 






0.38 ±0.12 


1.1 






1.78 ± 0.69 


J/"^ + Xc2 






0.69 ±0.13 


1.6 






0.44 ±0.17 


r]c + he 






0.308 ±0.017 








0.25 ±0.10 


he + XcO 






0.053 ±0.019 








0.075 ± 0.029 


he + Xcl 






0.258 ± 0.064 








0.132 ±0.051 


he + Xc2 






0.017 ±0.002 








0.004 ± 0.002 



pedagogical personnel of innovative Russia" (grant No. NK-20P/1). 

Appendix A: The coefficients Ai, Bi and Cj entering in the production amplitudes 

(16)-(18) 

These coefficients are the sums of the terms containing the parameters 
u = M^^j/{Mh^ + M^^j) and k = m/{Mh^ + M^^J. We present Ai, Bi and Q in 
numerical form using the observed meson masses and the mass of c-quark m = 1.55 GeV. 

e"^ + e- he + XcO 

Ao = 47.20 - 6.63r2 ± Coi(-4.66 ± ST.TOr^ - 8.05r^ ± O.Olr^) ± (Al) 
±Cio(-40.36 ± 55.54r2 - 6.32r^ ± O.Olr^) ± e.GSr^Cos ± C2o(3.24 ± 5mr^) + 
±^11(3.98 - 48.83r2 ± 53.37r^ - 7.63r^ ± 0.02r^). 

e+ ± e~ -> /ic ± Xcl 

Bi = -4r2 ± 2.73r^ ± Cio{3A2r^ - 14.22r^ ± 3.68r6) ± C2o{^A\r'^ - 2.73r^) ± (A2) 

±Coi(-2.75r2 - 11.24r^ ± 3.89r^) ± Cixi^l.'i'br'^ ± 25.47r^ - 26.77r^ ± 5.06r®) - 2.73r^Co2, 
^2 = -7.50r2 ± 2.73r^ ± Cio(24.13r2 - 19.67r^ ± 3.68r^) ± C2o(7.50r2 - 2.73r^) ± (A3) 
±Coi(16.45r2-17.04r^±3.89r^)±Co2(4.51r2-2.73r^)±Cii(-44.63r2+65.23r^-34.83r6±5.06rS), 
Bj, = 18.99 - 9.19r2 ± Cio(-51.66 ± 58.38r2 - 13.28r^) ± C2o(-18.99 ± 8.75r2) + (A4) 
+Coi(-28.91 ± 49.12r2 - 13.49^ ± Co2(-9.03 ± 8.21r2) + 
±Cii(85.84 - 164.44r2 + 109.27r^ - 17.99r^), 
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^4 = 4 - 1.82r^ + Cio(-3.42 + 9.48r^ - 1.84r^) + C2o(-2 + 1.82r^) + (A5) 
+Coi(1.51 + 7.53r^ - 1.95r^) + 1.82r^Co2 + Cii(-1.29 - 17.07r^ + 13.40r^ - 2.03r^), 

= Coi(2.50-0.12r^+0.02r^)+C2o(-2.88+0.44r^)+Cii(-2.13+0.28r^-0.08r^-0.05r^), 

(A6) 

= Coi(-2.48 + O.OGr^ + 0.02r^) + O.Ur'^C2o + (:7ii(2.12 - 0.22r2 + 0.14/ - 0.05/). (A7) 
e+ + e- hc + Xc2 

Ci = -2.36/ + Cio(8.74/ - 3.28/) + 2.36/^20 + Coi{0.35r'^ + 10.22/ - 3.22/) + (A8) 



+Coi(1.16 + 0.46/ + 0.41/) + Cii(-1.16 - 1.50/ - 0.71/ + 0.66/), 

Cg = -1.21 + 1.58/ + Cio(1.67 - 6.60/ + 1.64/) + C2o(0.42 - 1.58/) + Co2(0.82 - 1.58/) + 

(AlO) 

+Coi(0.46 - 6.65/ + 1.61/) + Cii(-0.33 + 19.72/ - 11.80/ + 1.73/), 
C4 = 1-21 + Cio(-1.67 + 0.78r2) - O.42C20 + Cii(0.33 + 1.70/ - O.6O/) + (All) 

+Coi(-0.46 - 0.38/) - 0.82Co2, 



Cs = Coi(0.47/ + 0.40/) + Cii(-0.5l/ - 0.57/ + 0.73/) - O.38/C20. (A15) 



Appendix B: The functions F-^ (r^) {k = 0, 1,2) entering in the production cross 



+Cii(-0.42r^ - 31.78r^ + 24.24r^ - 4.33r^) + 2.36/Co2 
C2 = -1.21 + Cio(1.67 - 0.78/) + O.42C20 + 0.82Co2 + 



(A9) 



C5 = Cii(-0.58/ + 0.89/) + O.8I/C01, 
Ce = -0.12/Coi + Cii(-O.Ol/ - 0.18/), 
C7 = Coi(-0.35/ + 0.33/) + Cii(0.42/ - 0.35/ + 0.55/), 



(A12) 
(A13) 
(A14) 



section (19) 



e+ + e -^hc + Xco 



Fo^°^ = 2.25 - 2.88/ + 0.68/ - 0.04r^ 



(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 



Ff' = -3.84 + 9.67r2 - 7.17/ + 1.43/ - 0.09/, 



Fr = 0.31 + 0.12/ - 0.50/ + 0.07/, 



^ = 1.64 - 6.16/ + 8.14/ - 4.33/ + 0.75/, 
fP = -0.44 + 4.09/ - 4.92/ + 1.38/ - O.ll/, 



^5^°) = 0.63/ - 0.72/ + 0.09/, 



F^"> = 0.02 - 0.38/ + 1.86/ - 2.12r'^ + 0.68r^ 
F^^^ = 0.76 - 9.05/ + 18.96/ - 13.49/ + 3.03/. 



e"*" + e -^hc + Xci 
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^({1) = 0.19 + 2mr^ - 4.32r^ + 2Mr^ - O.Slr^ (B9) 

^ = -0.32 - 11.48r2 + 30.85r^ - 28.58r^ + ll.OOr^ - 1.47r^°, (BIO) 

^ = -0.37 - 3.96r2 + 8.45r^ - S.lOr^ + 0.98r^ (Bll) 

Fg^^^ = 0.14 + IS.Tlr^ - 51.51r^ + 64.29r^ - 37.91r^ + 10.34ri°, (B12) 

^ = 0.14 - 6.78r2 + 21.01r^ - 22.91r^ + 10.03r*^ - 1.49ri°, (B13) 

= -2.13r2 + 5.02r^ - 3.81r^ + 0.92r^ (B14) 

Fg^^^ = 0.03 + 4.78r2 - 21.19r^ + 35.38r^ - 26.94r* + 9.03r^°, (B15) 

F^^i) = -0.24 + 38.00r2 - 133.92r^ + 186.34r^ - 123.58r* + 37.52r^°. (B16) 
e"^ + e~ hc + Xc2 

Fq^^^ = 0.23 - 0.81r2 + 0.98r^ - 0.40r^ (B17) 

^ = -1.72 + 6.75r2 - 9.73r^ + 5.83r^ - 1.13r^ (B18) 

F^^'^ = -0.47 + l.eir^ - 1.91r^ + 0.76r^ (B19) 

= 3.18 - 13.83r^ + 23.15r^ - 18.01r^ + 6.32r^ (B20) 

Ff ) = -1.92 + 7.56r2 - 10.88r^ + 6.34r^ - l.llr^ (B21) 

Ff ^ = -0.47 + 1.63r2 - 1.96r^ + 0.80r^ (B22) 

F^ = 4.14 - 17.68r2 + 29.23r^ - 21.88r^ + 6.98r^ (B23) 

Fi^^ = 13.29 - 58.23r2 + 97.89r^ - 75.59r^ + 25.74rl (B24) 

Appendix C: Effective relativistic Hamiltonian 

For the calculation of the relativistic corrections in the bound state wave functions \E'^ we 
consider the Breit potential (23). It contains a number of terms which should be transformed 



in order to use the program of numerical solution of the Schrodinger equation |39l| . The 



rationalization of the kinetic energy operator can be done in the following form [4& 



T = 2vW^=2^!±£,«l^ + ?|!. (CI) 
VP + ^ E 



where m is the effective mass of heavy quarks, 



m = ^ = )LLllL . (C2) 

vlff should be considered as a new parameter which effectively accounts for relativistic 
corrections in (CI). We take numerical value of vlff = 0-54 GeV^ for P-wave charmonium 
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TABLE III: The parameters of the effective relativistic Hamiltonian and masses of V-w&ve char- 
monium states. 



Meson (cc) 




Plff, GeV2 


m, GeV 


M*^, GeV 


j^exp^ GeV [42] 


XcO 




0.54 


0.857 


3.418 


3.415 


Xcl 




0.54 


0.857 


3.493 


3.511 


Xc2 


i^P2 


0.54 


0.857 


3.557 


3.556 


he 




0.54 


0.857 


3.499 


3.525 



states (see Table III). The second term in the Breit potential (23), which also has to be 
transformed, takes the form: 



AU 



3m^r 



(C3) 



It has the similar structure as the operator of effective kinetic energy from the Hamiltonian 



39 in order to include 



Hq. So, we change slightly the code of the Mathematica programm in 
the correction AU directly in the initial Hamiltonian. 

At last, there is a need to transform the spin-spin and spin-orbit interactions in (25) which 
have the 1/r^ behavior at small r. For the purpose of solving the Schrodinger equation we 
consider the regularization of such terms due to the account of the relative motion of heavy 



37| . The nonsingular potentials in both 
fi)/r^ with /i = (1 + 2mr)e~^'^'' and 



quarks which was discussed many times in |32l . l36l 
cases have the following structure at small r: (1 — ji)/ 1 wini ji 

/2 = (1 + 2mr + 4m^r^/3)e~^'"'' for spin-orbit and spin-spin interactions correspondingly. In 
Table HI we present the results of the calculation of the "P-wave charmonium mass spectrum 
and a comparison with the existing experimental data. The obtained masses agree with the 
experimental ones within an accuracy 1%. So we can suppose that the obtained effective 
Hamiltonian allows to account relativistic corrections in the bound state wave functions with 
sufficiently good accuracy. 
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